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Abstract 

We revisit the problem of decay of a metastable vacuum induced by the presence 
of a particle. For the bosons of the 'master field' the problem is solved in any number 
of dimensions in terms of the spontaneous decay rate of the false vacuum, while for a 
fermion we find a closed expression for the decay rate in (1+1) dimensions. It is shown 
that in the (1+1) dimensional case an infrared problem of one-loop correction to the 
decay rate of a boson is resolved due to a cancellation between soft modes of the field. 
We also find the boson decay rate in the 'sine-Gordon staircase' model in the limits of 
strong and weak coupling. 



1 Introduction 



The decay of a metastable vacuum state is a quite universal problem in quantum field theory. 
The decay proceeds through nucleation and subsequent classical expansion of the bubbles 
of the true vacuum. The classical bubbles can exist only starting from a certain critical 
radius at which the energy loss due to the surface terms is compensated by the gain in the 
volume energy. The formation of the critical bubbles is thus a quantum tunneling process P 
which tunneling can be described by an Euclidean-space configuration of the field, called a 
'bounce' [2]. The space-time nucleation rate of the critical bubbles, Wq, i.e. the probability 
of such nucleation per unit time and per unit volume, is proportional to the exponent of 
the classical action on the bounce configuration: wq oc exp{—Sci), while the pre-exponential 
factor requires a calculation of the functional determinant at the bounce. The exponential 
factor is readily found^l |2] in the so called thin wall limit, namely when the radius of the 
critical bubble is much bigger than the effective thickness of its wall. This limit is always 
realized at a small difference e between the vacuum energy density of the false and the true 
vacua, with the other parameter determining Sd being the surface tension fi of the bubble 
wall, i.e. of the boundary between the metastable and the stable phases. The pre-exponential 
factor in the bubble nucleation rate is known in a closed form only in (1+1) dimensional 
models j31|3], where in the thin wall limit it is determined only by the parameter e, with only 
partial results found in (2+1) dimensions [HI HE], and virtually no result known in the (3+1) 
dimensional case. 

Similarly to the behavior in the decay of a metastable phase in a thermal setting, the 
presence of matter in the false vacuum generally provides 'centers of nucleation' for the 
bubbles of the true vacuum. Thus one can consider the false vacuum decay induced by the 
presence of a particle OUni HI], by particle colhsions fTIH ITT]. by matter with finite densityjT^, 
as well as by the matter being in a thermal equilibrium where the problem goes back to the 
more conventional thermodynamic setting^j. In this paper we revisit the calculation of 
the bubble nucleation rate associated with the presence of a particle in the false vacuum. 
The particle-induced nucleation can also be viewed as the decay of the particle (albeit in 
the process the initial 'vacuum' state also gets destroyed), whose rate T generically can be 
written in the form F = Kwq, where the constant K, which can be naturally called the 
'catalysis factor', is the main subject of our consideration. The catalysis is most efficient for 
the particles which have zero modes localized on the boundary between the false and the true 
vacua. The reason for this behavior is that in this case the energy corresponding to the mass 
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of the particle m in the initial state is fully transferred to the bubble degrees of freedom, 
since in the final state the particle ends up as a zero mode localized on the bubble wall. 
This effectively corresponds to the upward shift by m of the energy at which the tunneling 
takes place [TUj. and results in K being proportional to the exponential factor exp(2mr), 
where r is the (Euclidean) time on the tunneling trajectory. The exponential behavior due 
to the shift of the energy for the tunneling trajectory can be found explicitly both in (1+1) 
dimensions|l()j and in higher-dimensional modelsf^. However the pre-exponential factor 
has been calculated only for the bosons of the master field in a (3+1) dimensional case[9], 
for which bosons the existence of the zero mode is always true. Here we calculate the pre- 
exponential behavior of the catalysis factor for the same bosons in lower dimensions, and 
also find a closed formula for this factor in (1+1) dimensions for a fermion, whose field has 
a zero mode on the inter- vacua boundary. The existence of such fermionic mode is a generic 
phenomenon and is guaranteed in the case where the mass term for the fermions changes 
sign across the bubble wallfHj. 

Our consideration, similarly to Ref.|U], is generally limited to models with weak coupling, 
which implies that the masses m of the both types of considered particles are small in 
comparison with the scale of the surface tension /x. In this case the deformation of the 
tunneling trajectory due to the energy shift by m^TIH ITT] can be neglected, so that, in 
particular, the tunneling time r coincides with the radius R (x fi/e of the critical bubble, 
r = _R, as it does in the spontaneous vacuum decay Furthermore, we also assume the 
applicability of the thin wall limit, which implies the condition mR^ 1, and which is always 
valid in the limit of small e. 

The catalysis factor K, as defined, has the dimension of the spatial volume. Thus it 
would be natural to compare the pre-exponential factor in K with the spatial volume of 
the critical bubble of the radius R. Under our assumptions we find that for a fermion in 
(1+1) dimensions this factor is indeed of order R, while the catalysis factor for the bosons 
is enhanced in comparison with the volume of the bubble by inverse powers of the (small) 
coupling constant. 

It should be noted that technically the bosonic catalysis factor in lower dimensions can 
be found by a straightforward application of the treatment of Ref.jHI. Such application is 
fully justified in a (2+1) dimensional case. However in (1+1) dimensions there is a potential 
complication in estimating the effect of the quantum fluctuations arising from an infrared 
behavior of the modes of the bosonic field over the bounce background. We demonstrate for 



2 



this case that the large infrared terms in fact cancel due to the specific properties of the soft 
modes. 

The material in the rest of the paper is organized as follows. In Sec. 2 we briefly review 
the calculation of the spontaneous decay rate of false vacuum and present a calculation of 
the decay rate induced by a boson of the scalar field, which defines the vacuum states. In 
Sec. 3 the problem of the infrared behavior of the one- loop correction to the calculated decay 
rate in (1+1) dimensions is considered and it is shown that this problem is resolved due 
to a cancellation of the contributions to this correction between the negative mode and the 
sum over the positive soft modes of the field of the bounce. In Sec. 4 the catalysis factor is 
calculated for a fermion in a (1+1) model. We then discuss the decay of metastable states in 
the sine-Gordon model with added linear term, the so called 'sine-Gordon staircase'. Using 
the equivalence of this bosonic model and the massive fermionic Thirring model in an 
external electric field, we find the induced decay rate for both the weak coupling limit (Sec. 5) 
and for the strong coupling limit (Sec. 6), the latter corresponding to a weak coupling in the 
Thirring model. Finally, in Sec. 7 we discuss possible implications of our calculation for other 
models. 

2 Spontaneous and induced decay of false vacuum 

In what follows we assume a situation where the energy density of a scalar field 0, the 'master 
field', has a local minimum at = 0+, which is higher than in a neighboring minimum at 
(j) = The vacuum state defined by the former minimum is referred to as the false vacuum, 
while the latter is the true vacuum. A typical example of such situation is provided by the 
well known model of a scalar field with the potential 



where A, and a are constants. At a = the potential has two degenerate minima at 
0± = ±f , while at small positive a the degeneracy is lifted in such a way that the minimum 
at 0+ has energy density bigger than that of by the amount e 2av. The vacuum state 
at 0+, being stable at a < becomes metastable at positive a and decays by nucleation and 
subsequent expansion of the bubbles filled with the phase 0_. At small a the surface density 
of the bubble wall can be approximated [H |2] by the surface density of the soliton with the 
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(2) 



interpolating between the two degenerate vacua in the hmit a ^ 0: 



J —V O 



(3) 



The mass of the scalar particles of the field propagating in either of the vacua is given (also 
in the limit a 0) a.s m = Xv. In a model with the total space-time dimensions equal to 
d the ratio m'^~^/fi coincides with the dimensionless coupling constant for the perturbation 
theory in this model. We assume throughout this paper that this ratio is a small parameter, 
which thus corresponds to weak coupling. 

In the Euclidean-space formulation of the problem of the false vacuum decay the cal- 
culation of the spontaneous decay rate amounts to a semiclassical evaluation of the imaginary 
part of the energy of the false vacuum from the path integral 



where the dots stand for other possible fields present in a specific model, Af is the normaliza- 
tion factor, and the integration is performed with the condition that the field approaches 
its false vacuum value 0+ at the boundaries of the space-time normalization box. The de- 
cay rate is then given by wo = 21m{lnZ)/VT, where VT is the space-time volume of the 
normalization box. 

The action functional S has a semiclassical saddle point at the configuration described by 
the bounce In the thin wall limit the bounce is an 0{d) symmetric bubble with the field 
0_ inside and 0+ outside, and the bubble wall, separating the two phases has the surface 
tension fi. The action for the bounce in this approximation is given by 



where Vb is the d dimensional volume of the bounce and Ab is its {d — 1) dimensional 
surface area. The action (0) reaches its extremum on a spherical bounce with the radius 
R = {d—1) fi/e, which is also the radius of the critical bubbles capable of classical expansion 
in the Minkowski space-time. 

The spectrum of small deformations of the bounce around the extremum contains exactly 
one negative mode, corresponding to an overall variation of the radius. This mode in fact 




(4) 



S = ij,Ab - eVB , 



(5) 
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gives rise to the imaginary part|2j of the path integral in eq.(j3]). Furthermore this spectrum 
also contains d translational zero modes, the integration over which introduces the factor of 
the space-time volume VT in the contribution of the bounce to the energy of the vacuum 
state. 

The decay rate of a particle of the field in the false vacuum can be calculated jHl El 
by considering the imaginary part of the contribution of a bounce to the Euclidean-space 
propagator of the excitations cr(x) = — 0+ of the field 0: 

D{x, y) = ^ J (r{x)a{y) e-^l-^'-l D0 . . . (6) 

in the limit of large separation L = \x — y\. Indeed, the contribution of the bounce to the 
correlator ©, as shown in Fig. la, has the generic form 

SD{x, y)=^-woJ d'^z F{x - z,y - z) Do{x - z) Do{y - z) , (7) 

where Dq[x) is the free-particle propagator in the vacuum 0_|_, satisfying the equation 

[-3'' ^ m") D{x) = b^'^\x) (8) 

and the factor («/2) w^d'^z is the proper measure of integration over the coordinate z of the 
center of the bounce, as follows from the consideration of the bounce contribution to (the 
imaginary part of) the vacuum energy. 

Let us consider the contribution to the integral (|7|) arising from the configurations, where 
the bounce is far (in units of its radius) from either of the points x and i.e. where 
|x — z| ^ R and |y — -z] ^ R. The propagators Dq{x — z) and -Do(y ~ z) in the integral in 
eq.© describe the exponential attenuation of the correlation {D{x) ~ exp(— m|a;|) at large 
separations, while the form factor F{x — z,y — z) does not have this exponential behavior. 
For this reason at \x — y\ = L ^ R the integrand in eq.© is maximized for z lying on 
the straight line running between x and y: z^ = s [xy — y^)/ and the integration can be 
split into the longitudinal, over the parameter s along this line, and the transversal, over zj^. 
The integration over z± can be done by the saddle point method, so that the form factor 
F{x — z,y — z) can be replaced by its value at z± = 0, and the essential configuration to 
be considered is the one shown in Fig. lb. As will be discussed few lines below, when the 
bounce is far from the endpoints of integration over s, i.e. s ^ R and L — s ^ R, the value 
of the form factor in fact does not depend on s and is a constant Fq. Since the contribution 
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a) b) 

Figure 1: The configurations for tlie bounce contribution to tlie propagator (eq.Q). A 
generic configuration (a) and the ahgnment of the bounce position (6), dominating the 
integral in eq.© at large \x — y\. 



of the excluded regions around the endpoints is only of relative order R/L, the integral in 
eq.(I7)) can be replaced at large L by 

5D{x, y)='-woFo J d^z D^{x - z) D^{y - z) , (9) 

where Fq should be calculated from the configuration shown in Fig. lb. 

The expression ^ for the modification 5D of the propagator by the bounce can be 
compared with the first-order correction to the propagator due to a small shift of mass by 
5m?, m? — i> m? + 5m?, in eq.®- In the standard way one finds 

5mD{x, y) = -5m^ J d^z D^{x - z) D^{y - z) . (10) 

Thus the contribution Q of the bounce to the propagator of the boson in the false vacuum is 
equivalent to an imaginary shift of the boson mass: 5m'^ = —{i/2)woFo, which corresponds 
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to the particle decay rate given by F = FoWo/(2m), so that the catalysis factor K is found 
as 

K = ^. (11) 

The factor Fq can be readily found [U^ for the discussed here case of the bosons of the 
classical field of the bounce. Indeed, consider the classical field just outside the bounce, 
i.e. at the distance r > R from the center, such that r — R ^ m~^, but still r — R <^ R. 
The former condition ensures that the field is described by its asymptotic approach to the 
vacuum value 0+, while the latter implies that in this region the curvature of the bounce 
wall can be neglected in a calculation of this asymptotic behavior. Thus one can consider 
instead the asymptotic behavior of the field in the limit e — > 0, i.e. of the field of the stable 
soliton separating two degenerate vacua. This asymptotic behavior has the form (j){x) — (f)^ = 
— 2f exp[— m (r — R)], where in the model described by the potential (0) v coincides with the 
corresponding parameter in the potential, while in a generic model v ~ (0+ — 0_)/2. On the 
other hand in the 0((i)- symmetric problem the asymptotic approach of the scalar field to its 
vacuum value is described by the solution of the linearized spherically-symmetric equation, 
equivalent to the homogeneous part of eq.®, and reads as, 

0(r)-0+ = CDo(r) (12) 

where the free boson propagator in d dimensions has the well known expression in terms of 
the modified Bessel function Ky{mr): 

^o(^) = (27^)^/2 r^/2-i -^'^/2-^(^^) • (13) 



The constant C in the asymptotic expression (fT^ is found by comparing the two expres- 
sion for 0(x) — 0+ in the discussed region just outside the bounce and using the standard 
asymptotic formula for the function Ky{mr). In this way one finds 

C= -4(27r)'^/2-l^(3-d)/2^(d-l)/2^gr„K _ (^4) 

Using then the expression (fT^ for the field with thus determined constant C, one finds the 
product of the classical fields a{x)a{y) in the integral in eq.® in the configuration shown 
in Fig. lb, corresponding to the constant Fq in eq.© given by 

Fo = = 16 (271)^^- W-'^ i?^'^-!) e^"^^ , (15) 
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which indeed does not depend on the position of the bubble along the straight line connecting 
the points x and y as long as both these points are sufficiently outside the bounce. The 
catalysis factor thus can be found from the relations and ()15|) in the form 

where Vd-i = n'^'^^^^^'^ R^^^ /T[(d + l)/2] is the spatial (rf — 1 dimensional) volume of the 
critical bubble. As discussed in the Introduction, it is natural to compare the catalysis 
factor with this volume. The result in eq. ()16j) shows that besides the classical exponential 
factor the catalysis is additionally enhanced by the factor m^~'^f ^ in the pre-exponent, which 
is the inverse of the small dimensionless coupling in the theory. 

3 Boson-induced decay in (1 + 1) dimensions 

The formula for the catalysis factor in eq.()16p reduces in (3+1) dimensions to the result of 
Ref.jO], and in other dimensions it presents a rather straightforward generalization. There 
is however one point, of a special importance to a (1+1) dimensional case, related to the 
effect of the quantum fluctuations on the essentially classical result in eq. ()l(i|l . Generally, 
the effect of the quantum fluctuations (the loop correction) is expected to be suppressed by 
a power of the coupling constant as compared to the classical contribution. In the discussed 
problem this expectation is true at d > 2, however in a (1+1) dimensional problem this 
expectation is potentially jeopardized by an infrared behavior. Indeed the eigenvalues of 
the second variation of the action for the fluctuations of the shape of the bounce, described 
by the effective action are proportional to R~^. The modes with these eigenvalues are 
localized on the bounce boundary and describe the soft part of the spectrum of the modes 
of the field around the stationary bounce configuration, as opposed to the modes, whose 
eigenvalues start at 0(m), and those 'hard' modes describe the excitations propagating in 
the bulk as well as possible deformations of the proffie of the field across the bounce wall. 
Let us estimate the contribution of an individual soft mode with the eigenvalue Cn/R"^ to the 
correlator (jHl), with c„ being a number. All such modes originate from local shifts of the wall 
of the bounce, so that the field proffie of an individual mode is proportional to the radial 
derivative of the field of the bounce, 0'(r). The field an of a normalized to one mode in (1+1) 
dimensions is then parametrically estimated at the distance r > R, such that r — R^ m^^, 
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but still r — i? ^ i?, as 



^^{^r) ~ ^ e-"^(^-^) ~ e-"^('^-^) , (17) 
Y /iix V R 



where it is taken into account that J{(j)')'^dr ^ fi mv'^, and any numerical factors are 
dropped for a parametrical estimate. The contribution of such mode to the correlator (jH)) is 
then proportional to 



— cr„lri)a„(r2) ~ c"^ (m R) e^"^^ e-"^(^^+'^^) 



In this estimate the large factor (m R) stands in place of the factor f ^ in the similar product 
for the classical field. Thus the magnitude of an individual mode contribution to the corre- 
lator relative to the classical part is described by the parameter which in-spite of the 
expected suppression by the dimensionless coupling f ~^ is infrared unstable at large R. 

We will show however that the sum over the soft modes gives zero for this infrared 
contribution due to a cancellation between one negative and all positive modes. In order to 
demonstrate this we consider the parametrization of the shape of the bounce in the polar 
coordinates (r, 6) on a plane, so that the effective action (0) for the soft modes takes the 
form 2 

S= rU^^V^:^T^-lerA d0 = ^+ i{p''-p'')d9 + 0ip''), (19) 
JO V 2 / e Jo 2 

where the latter expression shows the two first terms of expansion in the small deviation 
p = r — R of the radial variable r from its stationary value R = p/e, and the dot stands 
for the derivative over 6. The quadratic part in this expression has one negative eigenmode 
p = 1/ \^27i and the spectrum of zero and positive double degenerate eigenmodes: 

p^^^ = cosnO , and p^^^ = ^ sinra^ ; (ra = 1, 2, . . .) . (20) 



The spectrum of the eigenvalues is proportional^ to (n^ — 1) with the negative mode corre- 
sponding to n = 0. 

Let us consider now the configuration shown in Fig. lb with the bounce located on the 

line connecting the points x and y. Let the angle 6 be defined as measured counterclockwise 

from the downward vertical connecting the center of the bounce with the point x, so that 

^The proportionality coefRcient is not important for this discussion. It can be noted however that in 
terms of the normahzed modes for the field (j> the eigenvalues are (n^ — 1)/ E? . 
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6 = TT corresponds to the upward vertical connecting the same center with the point y. 
Clearly, the contribution of the fluctuations of p to the propagator is proportional to 

(Ho)+pwnc.E '""'°,'/''f"' . (21) 

n n I 

Note however that the sum p(0) + p(7r) is not vanishing only for the negative mode and for 
the positive modes of the first type, with even n, i.e. n = 2k. Thus the sum in eq.(j2H) 
is proportional to the numeric sum 

1 OO 1 

-2\S4F^-0. (22) 

k=l 

where the first term is due to the negative mode and the sum runs over the positive modes. 
The arithmetic identity (j22|) explicitly demonstrates that the infrared contribution in a (1+1) 
dimensional model cancels between the negative mode and the sum over the positive ones. 

Let us also remark on the decay of a moving particle in the false vacuum. If particle 
moves with a constant velocity then the probability depends on the velocity through the 
standard Lorentz factor. However if it moves with the constant acceleration situation is 
more subtle since in the particle frame vacuum behaves as the thermal bath due to Unruh 
effect. The effective temperature is defined through the acceleration as 

n,j = ^ (23) 

hence probability of the particle decay is modulated by the thermal effects. The most 
essential effect corresponds to the possible deformation of the classical bounce. Since tem- 
perature corresponds to the periodicity in the Euclidean time then the deformation of the 
bounce happens when period corresponding to the temperature becomes comparable to 2R. 
That is deformation of the bubble emerges if the acceleration is larger then 

vre 

acrit = —■ (24) 

and our approximation fails. 



4 Fermion- induced decay in (1 + 1) dimensions 

The very existence of fermions in a model is known to modify (without any fermions being 
present in the initial state) the pre-exponential factor in the rate of the false vacuum decay 
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in the situation where the complex fermion field ip has a zero mode on the boundary between 
the vacua (in the limit e — 0). Such situation takes place when the mass term for the fermion 
changes sign between the two vacua^j. In particular the rate wo for the spontaneous decay 
of the false vacuum in (1+1) dimensions receives a factor of 2 in comparison with purely 
bosonic theorv fTBl IT7]. This doubling corresponds to the existence of two final states in the 
false vacuum decay in (1+1) dimensions viewed as a spontaneous creation of a kink-antikink 
pair: one state where both the kink and the antikink are created with the fermion zero mode 
empty, and the other state is where is a zero-energy fermion on the kink and a zero-energy 
antifermion on the antikink. 

In what follows we assume that the interaction of the fermion field with the scalar field 
of the bounce is such that there exists a zero fermion mode on the kink separating the two 
vacua. In order to find the effect of a fermion on the probability of nucleation of a critical 
bubble we consider the bounce contribution to the fermion propagator G{x, y) = {4'{x)ip{y)) 
in the configuration shown in Fig. lb. Clearly, the exponentially enhanced factor exp(2 mj i?), 
where is the mass of the fermion, arises from the contribution of the zero mode, whose 
propagation from the lower point of the bounce to the upper one does not contain any 
exponential attenuation. Assuming for definiteness that the fermion mass is positive (and 
equals nif) in the false vacuum, and choosing the 7 matrices as 71 = ui and 72 = <72, one 
finds the solution for the Dirac equation for the field of the zero mode of ip in the background 
scalar field (t>(r) of the bounce, 



M, + m(</.)]Vo = , (25) 



in the form 



-expj-y^ m[0(rO]rfr'|x(^) ^^./^ ' (26) 



where x(^) is a one-dimensional fermion field living on the bounce boundary and (nominally) 
depending on the length parameter i = R9 along the boundary. Notice that the classical 
equation for x reads x = 0. Finally, the constant C/ in eq. ()26|l is the normalization factor 
relating the normalization of ip and x ^iid satisfying the condition 



2C| y" exp -2 m{4>)dr' 



dr = l . (27) 



In what follows we rather use a related factor Cf defined as the coefficient in the expression 
Cf exp |- m[0(r')] rfr'j ^ Cf exp [mf{R - r)] , (28) 
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which is vahd sufficiently far outside the bounce where rrimin (r — i?) ^ 1 with m^m being 
the minimal mass scale in the model. Generally the factor Cf can be estimated as 

where / is a dimensionless function of the ratio of to masses of other particles in the 
false vacuum. In the limit where mj is much smaller than other masses one has /(O) = 1. 
In the model, where the scalar 'master field' is described by the potential (P) and the mass 
of the fermion is proportional to 0, the function / can be found explicitly: 

_ 2- r(n + l/2) 

The contribution of the fermion zero mode on the bounce the configuration shown in 
Fig. lb can be written, using the asymptotic behavior of the zero mode ()26|) . in terms of the 
one-dimensional propagator of the field x otl the boundary, g{£i,£2) = (x(^i)x^(^2)), as 

5Gix, y) = -'-^dh Cj e'^f^'R , ^ (1 + ai) g{0, nR) . (31) 



\x — z\\y — z\ 

Notice that for a complex fermion there is only one path for propagation along of x along 
the boundary from the bottom of the bounce to its top (assumed here for definiteness to be 
counterclockwise in terms of Fig. lb), corresponding to the final state, where the fermion is 
a bound state localized on the kink. The other path (clockwise) would be relevant for the 
vacuum decay induced by an anti-fermion, which in the final state is localized on the antikink. 
The expression in eg. 1)311) contains an extra factor 1/2, due to the fact that, as mentioned 
before, the spontaneous nucleation rate Wq in the theory with fermions contains extra factor 
of 2, due to the existence of two final states in the decay, which is to be compensated 
in the proper measure of integration over the coordinate of the center of the bounce d'^z. 
The propagator g has a very simple explicit form in terms of the sign function: (7 (£1,^2) = 
(1/2) sign(£i - £2), so that ^(0, nR) = -1/2. 

The expression in eq. (j!?T| can now be compared with the corresponding change of the 
free propagator Go under a shift 6m f of the fermion mass: 

lyjj^ g^l^' y\ 

SmG{x, y) = —Snif Sz Go{x — z) Go{z — y)—^ —Smrd^z — (1 + o"i) , (32) 

47]- I' 



\x — z\\y — z\ 

where the asymptotic expression takes into account the explicit form of the free propagator: 

Gq{x, y) = 77- {-Oidi + m) KQ{mj \x - y\) . (33) 
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Using this comparison and eq.()29|) one finds the imaginary part of the fermion mass shift 
corresponding to the decay rate of the fermion 

r/ = ^/(^) /?^oexp(2m,i?) = ^/(^) exp(^-^ + 2m,i?j . (34) 

Here in the latter transition are used the exphcit expressions: wo = (c/tt) exp(— Tr/i^/e) and 
R = fi/e. One can readily see that in the fermion expected, the pre-exponent in the 

catalysis factor, Kf = {tt/2) f{mf/m) R exp(2m/i?), is indeed of the order of the spatial 
size of the critical bubble. 



5 Meson decay in weakly coupled sine-Gordon model 

In order illustrate the universality of the derived results let us discuss one more example 
where the particle decay in the false vacuum happens in two dimensions. We shall derive the 
probability of the decay of the electrically neutral meson bound state in the Thirring model 
with the pre-exponential accuracy. 

Consider the sine-Gordon theory with the Lagrangian 

1 ry 

LsG = ^{d<pf + j^cosm (35) 

and add term (e/3/27r) which yields the situation with the metastable states. This theory 
upon the two-dimensional fermionization is equivalent to the massive Thirring model with 
the Lagrangian 

Lth = i^d^Yi^ - ^9 fju + /i + ^0 Jo (36) 

where £^ = (1 + f )^\ and = ip'ju'ip- One can identify /i with the sohton mass in the 
sine-Gordon model, and d^A^ = e. In what follows we shall assume that < Air which is 
the condition for the bound state of fermions to exist in the Thirring model. The solitons in 
the sine-Gordon model get mapped into the fermions in the Thirring model while the field 
gets mapped into the fermion-antifermion meson bound state. 

The linear perturbation term in the sine-Gordon model corresponds to the constant 
electric field in the Thirring model realization, so that the problem of false vacuum decay 
can be discussed in both formulations. In the Thirring model it corresponds to the Schwinger 
pair production. The probability of the spontaneous vacuum decay in the sine-Gordon model 
and equivalent Schwinger process in Thirring model has been found in The one-loop 
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result coincides with the general formula wxh = (e/27r) exp(— vr/i^/e), while in the special 
case of = in, corresponding to g = 0, the exact result is found as 

WTh = ln(l - e--^'/^) (37) 

Note that in this case the bose-fermi equivalence allows to perform the summation over the 
multiple bounces in the sin-Gordon theory. 

Now we can discuss the decay of the false vacuum in the presence of a particle correspond- 
ing to the field in the sine-Gordon model. In the weak coupling regime for the bosons, 
i.e. at small (3, the soliton is much heavier then the boson particle, which corresponds to 
the situation where the external particle does not deform the classical bounce configuration. 
Hence the decay rate can be immediately read off eq. (|T6|) . In this case the process corre- 
sponds in the Thirring model to the nonperturbative decay of the light electrically neutral 
meson in the electric field and the catalysis factor of this process is 

KTH = f,^e'-''^/^ (38) 

where mi, is the meson mass. Note that this process is the two-dimensional counterpart of 
the induced Schwinger processes discussed in four dimensions in ^HJ in the exponential 
approximation. 



6 Meson decay in strongly coupled sine-Gordon model 

The boson-fermion correspondence in this model actually allows to find the meson decay rate 
in the limit, opposite to what has been considered so far in this paper, namely for strongly 
coupled bosons, when the boson mass is close to the kink-antikink threshold. This limit 
corresponds to a small positive g, and the boson mass (at e ^ 0) is = 2 fi — fig^. The 
near-threshold dynamics of the soliton-antisoliton pair can be considered nonrelativistically 
as a motion of a pair with the reduced mass /i/2 in the local potential U{x) = —2g6{x), 
which correctly reproduces the energy of the bound state (the boson). For a nonzero e the 
nonrelativistic Hamiltonian for this system takes the form: 

H = — -ex-2g5{x) . (39) 

The problem of the boson decay in the false vacuum is reduced in terms of this equivalent 
nonrelativistic system to that of ionization of the bound state in the external electric field e. 
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In order to solve the ionization problem we start with considering the Euclidean time 
propagator ("the heat kernel") defined as /C(x, r) = {x\ exp(— if r)|?/), and the correspond- 
ing energy dependent Green's function at the negative (i.e. below the threshold) energy 
E = -t^/^i: 



G l^x, y; -—j = IC{x, y; r) exp rj dr (40) 

at X = and y = 0. We remind that if only the kinetic term is retained in the Hamil- 
tonian (jSHl), i-e. at e = and g = 0, these functions are /Co(0,0;r) = (47r r//i)~^/^ and 
Go(0,0; —ti^/fi) = fi/{2K). At e = and a nonzero g the Green's function for the Hamilto- 
nian (fH^ is found as 

r rnn ^^/,\ ^0(0, 0; -kV/x) k 1 .... 
G,=o(0, 0; -/€ //i) = . 2TT = ^1 r- 

The latter expression contains explicitly the pole at k = ^g corresponding to the bound 
state. 

When both the e and g are nonzero the equation 1)4111 is modified by replacing the Green's 
function G'e=o by that for a nonvanishing e: G^. The latter Green's function can be expressed 
in terms of the corresponding propagator /Ce(0, 0; r) which can be found in the textbook 



exp 



4:71 T 



fi{x — yY e{x + y) 



r H 



4r 2 12 /X 



(42) 



a f 0, 0; = r exp i^r'-- r) dr , (43) 

and the pole position is thus determined from the equation 

2^G,(^0,0;-^j =1 . (44) 



The peculiarity of the latter equation is that the integral in eq.()43p is formally divergent, 
which is the usual situation in a calculation of the energy of an unstable state. In order to 
make physical sense, both that energy and the integral in eq. ()43p should be understood as 
a result of an analytical continuation in the parameters of the model from the region where 
the considered state is stable. In terms of eq. ()43)1 this corresponds to a continuation from the 
region of (formally) negative where the integral is convergent. The result of such analytical 
continuation to physical positive can be formulated as follows j2|: The integration runs 
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Figure 2: The contour of integration for the integral in eg . (1431) . 



along the real axis of r from r = to the value of r where the integrand has minimum, i.e. 
to tq = 2k,/ e. From that point the contour of integration should be turned parallel to the 
imaginary axis of r, corresponding to the direction of the steepest descent (see Fig. 2). This 
contour rotation gives rise to an imaginary part of the integral, and hence to an imaginary 
part of the energy of the resonant state, corresponding to the decay width of the resonance. 
Following this procedure one can readily find the real and the imaginary parts of the integral 
in eq. pHjl and reduce the equation (jl^ for the position of the pole to the form 



It can be noted that the exponential factor in this formula is the standard WKB tunneling 
exponent in a linear potential, while the pre-exponential factor is a new result. The described 
derivation of the formula assumes that the integral in eq. ()43p can be evaluated in the 
saddle point approximation, which implies that the parameter in the exponent in eq . (j46|) is 
large, i.e. that ji^ ^ e. 

7 Discussion 

In this paper we have refined the calculations of the decay rate of the boson in the false 
vacuum and have found the decay rate of fermion in the false vacuum in (1+1) dimension 




(45) 



which corresponds to the decay rate of the bound state 




(46) 
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with the pre-exponential factor. All calculations, except for the one in Sec. 6, have been 
performed in the approximation when the back reaction of the external particle on the 
Euclidean bounce solution can be neglected. 

The account of the back reaction amount to several new effects which are different for d=2 
and d > 2. In the (1+1) dimensional case the back reaction deforms the classical solution 
which deformation has been described classically in Ref.^U] however the calculation of the 
pre-exponential factor is beyond our approximation. Such calculation could be potentially 
interesting from the stringy perspective. Indeed, the worldsheet theory on nonabelian string 
in several models can be identified with the CP^ model (see for a recent review) which 
has one true vacuum and a set of metastable ones. At large N this theory can be treated 
perturbatively and the issue of the decay of metastable vacua or in other terns exited strings 
can be discussed. Similarly one can discuss the fate the different excitations on the exited 
metastable string which is just the problem we have considered. In some situations the 
pre-exponential factor is of the prime importance since in some range of parameters the N 
dependence disappear from the exponent |2Tj. However it is unclear if the regime with the 
negligible back reaction could exist in the worldsheet theory. It seems that the analysis 
similar to one in Section 6 could be applicable in this case. 

In higher dimensions the situation is more complicated. The point is that the initial 
particle evolves along a trajectory in the complexified Minkowski space^H]. At the first 
stage of the process the initial particle "produces" the oscillating bubble in the Minkowski 
space which later develops the path in the Euclidean space. The overlap of the initial particle 
and the bubble happens out of the real axis. Let us also note that there is a possibility of 
the resonant decay of the particle in the false vacuum when the particle mass coincides with 
the energy levels of the quantized bubbles in the Minkowski space. 
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